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system ano 1 ngui r i no as to the nacur e v_>l i ts 1 1 mix t sets 
C at.tr act or st in various cases Cin the spirit of CS 33. 

Thie iteration Cl. 13 is considered in a different context as 
a special case of the CR n -»[R n map 



V = A*V + B*V /IIV II Cl. 33 

n+1 n n n 

in C101. based on work in [2] Calso reported in C 333. Further 
deta i 1 s on the i ter at i on C 1 . 33 may be f ound in [10,12,133. and 
the forthcoming [43. Although C 1 . 33 is only a linear 

perturbation of the well-behaved iteration Cl. 13, it exhibits 
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The g o vp p m^r* t h c^d C 1. . cC' 1 h cas the pr **jpe r t y t hei. t i 2 A is ca 

nondefective matrix with a dominant eiaenva.lue, say X , and ^> 

i o 

has a rionzero projection ori an 01 oenvactor associated wi th 
this dor, 1 riant eigenvalue, then 

Li — > A. as n — ^ CD 

n 1 

and v converges to an eigenvector associated with X and with 

n ~ 1 

unit t norm. If v does not have a nonzero component along 
a> o 

an eign n ve cior associated with X and infinite precision 
an thine - ic is used then we must consider the eigenvalue of 
largest modulus along which does have a nonzero component. 

In actual computations, however , a component along an 
eigenvector associated with X^ would almost certainly be 
introduced eventually and magnified in successive iterations 
[1, p. 1 45] . The same results are found if A is defective 

([considering now' principal vectors [7, p. 3] rather than just 
eigenvectors!) but the convergence is much slower. 

When the dominant eigenvalue is real, the method converges 
to a fixed point. However, when the dominant eigenvalue is a 
complex conjugate pair, the method generally fails to 
converge. Methods exist to recover information in such cases 



[5, p.£o7] but they tend to be somewhat involved. 
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3. Conic Orbits 



If t is nonsinaular then the points V > V > V . . . . of the 

1 2 3 

orbit of Cl.il all lie on a conic def i ned by the matrix 

G = C B*B T D 1 
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for every n>l and for any V which is nonzero. Clearly G is 
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positive definite symmetric, and so the points V , V , . 

must all lie on the hyperellipse defined by 



V *G*V = 1 
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If B is singular, a similar result holds. In order to 

handl e si mu 1 taneousi y both the case wh e r e B is si mpl e and the 
case where B is defective we state the result in terms of 
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Theorem 1: Suppose B is a real square mxm matrix with 0<q<m 

nul 1 eigenvalues. If q> 0 and has a nonzero component along 

a principal vector associated with a nonzero eigenvalue, then 
all orb ts of the power iteration Cl. ID are constrained to a 
hyper ell ipse in C m-qD -di mensi ons C f or ail but finitely many 
nD . Otherwise, the orbit reaches the origin in finitely many 
.ti c ns . 

Proof: First, note that 

V = B*V /liV II 

1 o o 

V = R * « y ii V I 
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J = R *B*R 

be the Jordan normal form of B for some nonsingular R. 
Substituting this into C 3 . Ej gives 

V = CR*J*t? -1 :> ri *V /IIB r ' _1 *V II 
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= R*J n *R _1 *V ^ il B r ' _1 >tV II 
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a = l/IIB n **V II 
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is a scalar C f or each nD . Clearly, if Z has no nonzero 
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component along a principal vector of J Cequi valently, if 

has no nonzero component along a principal vector of B J that 

is associated with a non-null eigenvalue, the term J *Z in 
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undfei' 1 nv d and th^ i idralxon slops. Olhyrwi 9^, for n 

sufficiently large Cit suffices that ndriD , all Jordan blocks 
in J associated with a null eigenvalue will have become blocks 
of entirely zeros in J n Csince these Jordan blocks are 
nilpctentD. Now, the principal vectors belonging to a given 
Jordan block do no! interact with the remaining principal 
vectors, in the sense that if v is a principal vector 
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associated with J CX), a Jordan block of J, then J n *v 



* J 

involve! only a linear combination of principal vectors of J 

the. t ar» also associated with J CXk Therefore for n large 

enough ’hat ail nilpotent Jordan blocks have become entirely 

zero submatrices, the vector 
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can be written in terms of a basis consisting of only the 
rernainirg ( m-c) principal vectors. Thus the iteration lies in 
a C m-q") -di mens j onal subspace of CR™, and a suitable change of 
variables car then be used to transform the iteration into one 
of the form 



V = C*V I! V II 

n-*- i n n 

where V* is a C m-qj - vector for every n and C is a real 
nonsingular Cm-qDxCm-qD matrix. Hence, in this subspace, the 
iteration is constrained to the hyperellipse determined by the 
ma t r i x 

T -1 
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Cas was shown previously for the nonsingular case) for all but 
f 1 ni tel y many n . □ 



We emphasize that this is not an asymptotic result; after a 
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hyper ©1 ] 1 ps^ C assuming infinite precision) . We now wish to 
look at the orbits on the attract.! ng hyperellipses. 
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Suppose that B is nondefect ive and nonsingular, and let 
be a set of eiaenvectors associated with the 
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where a is the component of along x . Computing the power 



iteration Cl. 13 by C3.23 gives 
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^ l t t i ^ v l i 

1 = 1 J V = i 



C 4 . 1 J 



Of nil . 



Now suppose that B has a (repeated!) dominant eigenvalue a 
with multiplicity l<r<m, that is, 

x = x = • • • = X 

12 r 

and 



IX | > |X | > 
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> IX 



Removi no a factor of X n from the numerator and X° 1 from the 

i i 



denominator of C4.13 gives 

V = X n - [a X + 

n 1^11 
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X * j a x + 
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r = x -ex y\x p n 1 C4. 

r> 111 

= C sgnC X j r ' 4 * X 

^ i i 

C provided that at least one of c* , . . . , c* is nonzero; 

i i 

otherwise, we begin the analysis anew bv considennc X 0. 

r + 1 

Hence ii X >0 the iteration tends to a fixed point, and if 

i 

X^<0 the iteration tends to a symmetric Cin the origin} period 
two cycle on the points ±X^x. 

Now suppose that B has a real positive eigenvalue X^ of 
multiplicity p, so that 



X = X = 

i ? 



= A 



and a real negative eigenvalue -X^ of multiplicity q, so that 

= X 

p + 2 p -*-0 

and suppose further that 

■ • • ^ I A. | 
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| X I / j X 
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f F p * q r °l 
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r, 1 ^ \ i ^ \ \ p + q^i 1 J 
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J Vy'hte'T tr the — IS puS 1 ti ve when n IS 

r ) 

even and negative when n is odd in the numerator, and 

contrari ly in the denominator. In the limit, we have that 

C approx: matelyj 

v = r *cy + y'J/Wy - y II 

n n 1 2 l 2 

V = v -Cv - y D/lly + y II 

n+1 ‘ n+1 "1 1 2 
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1 V L 
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y = £ & x 
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Since A is positive, y =X for all n, and so we have a period 

1 n 1 
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V = X ’Cv + y ! / II y -vl 

o 11 2 1 2 

V = X -Cv - v D/llv + y II 

rt +1 1 ^ 1 2 ~ 1 2 



which in general is not a symmetric orbit Cin the origin^. We 
will not in general have II V ll =X in this case, but note that 

n 1 



II V ll = |X | • II y 

n 11 



V II A M V “ V I 
^2 ~ 1 2 



II V II = |X I • If v - y II Ally + y I! 

n+l 1 ' 1 2 *1 ^2 



so that 



li V II • II V ii = IX I 

n r >+ 1 1 

in the limit of large n; the modulus of the 
ei genva.l .je is the geometr i c mean oi the norms ot two 
success! vp i terates Ci n the 1 i m: 1 1 . Note that the asymmetry 

of the orbit allows us to distinguish this case from the case 

a l ways gi ves 

_ r oi i.j iTtv^s t the ti me C x , tr . wh^n the 

esul ti r g asymptotic orbit is indeed asymmetric!. 

Now' J f t us suppose that E has a complex conjugate pair of 
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, . , f -v 

= A | o< e x + a e 

1^1 1 2 
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x + OCX /X ! 
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v<n— l >$ -v(n-l><9 N ^ n-1 
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where £-'argCX J . In the limit, this becomes 



V = 

n 



, . , f -in G 

= X • I a e x + a e x e 

1^1 1 2 2 J 

v(n— 1>© -Un-1>© 

o e x + a e x 

1 12 2 
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Thus if 9 i s such that expC \9~J =expC vC n+1 } 0! for some n Ci.e. 



if X^ and X^ are both in K for some n! then we have an 



the orb: t is aperiodic on the underlying hyrereli ipse. If the 
orbit is period n then we note that the geometric mean of the 
norms o: n consecutive iterates tends to | X |, as is easily 
seen by w. 1 ting out t ne proouct ol n iterates and noting the 
cancel 1 a t i on . 

In a similar way, if E has two pairs of complex conjugate 
eigenvalues of equal modulus and all other eigenvalues of E 
have lesser moduli, then we have a situation like that of 
C 4 . 3 f save that there is an additional angle to be considered. 
Henc< if the f j rsd pair alone would give a periodic orbit of 
period n, and the second pair alone ’would give a periodic 
orbit of period q, then the orbit of the iteration will be 
periodic with period n=lcmCp,qJ. If either pair alone would 
give an aperiodic orbit , then the orbit is aperiodic. Again 
the georevric mean of n consecutive iterates tends to | | 
when t h orbit is as vmpt oti cal 1 y period n. The obvious 

generalisation holds for more than two complex conjugate p^airs 
of equal moduli. In particular, in an even-di mensi onal space, 
say of dimension Er , choosing B to be a 2rx2r matrix with r 
complex conjugate pairs of eigenvalues, all of equal modulus 
and such) that 

x n e IF 

for ali i =1 , . . . , Er and n=l ,2,3, . . . gi ves a method for 
generating a sequence of points on the ellipse given by C3.1D 
by iterating C 1 . 1 J with an arbitrary nonzero V . 

If in addition to some number of complex conjugate pairs of 
eigenvalues with equal modulus there are some number of 



i n 



f-\ — « A , , 

s_x w -l C A- * ~ c 3 - 



.•— f , ,-W. | *— -. ,*«, — A * ^ ♦- ■ ■ f- ,•— ' - 1 1 * * >■ T ♦ ♦ ' — * Kw '- A A t*'. 

iivw* . liCMai-A > c i Cul CAUCllvai ^ i i *= 



A V. •-•. >_ w /-•. .•- , -4 



^ai,ic jiu 



. .*4 • ♦ . « • *— 

^ W- , 



then it is immediate that the iteration is aperiodic if any 
com pi ex conjugate eigenvalue gives an aperiodic orbit, and 
pe r i oci i c of her wi s e , wi t h per iod given by the least c ommon 
multiple of the individual complex conjugate pairs and the 
period t wo due to the negative real eigenvalues, if present. 

We have established the f ol 1 owi no theorem: 
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V — B*s \ r / ii \ I! 

1 n n 

where B is a real mxm nondefective matrix and V is a Given 

o 

nonzero m- vector . Suppose that V has a nonzero component 
along an eigenvector of B which is associated with an 

eigenvalue of B with maximum modulus. Then all iterates for 
n>rri are constrained to an ellipse in a subspace of B " of 
dimension rank C ED , and the asymptotic behaviour of the 
iterates is as follows: 

i J if B has a multiple real dominant eigenvalue, then the 

iteration tends to a fixed point; 

iij if B has a multiple negative real eigenvalue, then the 

iteration tends to a symmetric period two cycle; 

iiO if B has both positive and negative real dominant 

eigenvalues, then the iteration tends to a period two cycle; 

iv) if B has a complex conjugate dominant eigenvalue, then 
the iteration tends to a period n cycle if the eigenvalues are 
a multiple of an n th root of unity and is aperiodic otherwise; 

vD if B has multiple complex conjugate dominant eigenvalues 
and multiple real dominant eigenvalues, then the iteration 
tends to a periodic orbit with period equal to the least 



11 



